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We study the classical and absolute stability of Q-balls in scalar field theories with flat potentials 
arising in both gravity-mediated and gauge-mediated models. We show that the associated Q-matter 
formed in gravity-mediated potentials can be stable against decay into their own free particles as 
long as the coupling constant of the nonrenormalisable term is small, and that all of the possible 
three-dimensional Q-ball configurations are classically stable against linear fluctuations. Three- 
dimensional gauge-mediated Q-balls can be absolutely stable in the "thin- wall-like" limit, but are 
completely unstable in the "thick- wall" limit. 



I. INTRODUCTION 

Q-balls have recently attracted much attentions in cosmology [l| and astrophysics 0, [3, 0] ■ A Q-ball 
@ is a nontopological soliton jBl whose stability is ensured by the existence of a continuous global charge 
Q (for a review see 0, H, S flOl fT]| and references therein) , and a number of scalar field theory models 
have been proposed to support the existence of nonto pol ogical solitons. They include polynomial models 
@, Sine-Gordon models parabolic-type models |l3j |. confinement models [13, E, [l3| , two-field 
models d, [13] , and flat models [l| . 

From a phenomenological point of view, the most interesting examples are probably the supersymmetric 
Q-balls arising within the framework of the Minimal Supersymmetric Standard Model (MSSM) , which 
naturally contains a number of gauge invariant flat directions. Many of the flat directions can carry baryon 
(B) or/and lepton (L) number which is/are essential for Affleck-Dine (AD) baryogenesis [l^. Following 
the AD mechanism, a complex scalar (AD) field acquires a large field value during a period of cosmic 
infiation and tends to form a homogeneous condensate, the AD condensate. In the presence of a negative 
pressure (20l . [2]| , the condensate is unstable against spatial fluctuations so that it develops into nonlinear 
inhomogeneous lumps, namely Q-balls. The stationary properties and cosmological consequences of the 
Q-balls depend on how the Supersymmctry (SUSY) is broken in the hidden sector, transmitting to the 
observable sector through so-called messengers. In the gravity-mediated [23 | or gauge-mediated scenarios 
[H , the messengers correspond respectively either to supergravity flclds or to some heavy particles charged 
under the gauge group of the standard model. 

Q-balls can exist in scalar field potentials where SUSY is broken through effects in the supergravity 
hidden sector [2^. These type of Q-balls can be unstable to decay into baryons and the lightest super- 
symmetric particle dark matter, such as neutrafinos [131, gravitinos [1^, [2g| and axinos [131 ■ Recently, 
McDonald has argued that enhanced Q-ball decay in AD baryogenesis models can explain the observed 
positron and electron excesses detected by PAMELA. ATIC and PPB-BETS [l^. By imposing an upper 
bound on the reheating temperature of the Universe after inflation, this mode of decay through Q-balls 
has been used to explain why the observed baryonic {^b) and dark matter {^dm) energy densities are 
so similar i.e. noAi/^b = 5.65 ± 0.58 

Scalar field potentials arising through gauge- mediated SUSY breaking tend to be extremely fiat. 
Using one of the MSSM flat directions, namely the QdL direction (where Q and d correspond to squark 
fields and L to a slepton field), which has a nonzero value oi B — L and therefore does not spoil AD 
baryogenesis via the sphaleron processes that violate B-\-L [s^l , Shoemaker and Kusenko recently explored 
the minimum energy configuration for baryo-lcptonic Q-balls, whose scalar field consists of both squarks 
and sleptons [11] • It had been assumed to that point that the lowest energy state of the scalar field 
corresponds to being exactly the fiat direction; however in [11] , the authors showed that the lowest energy 
state lies slightly away from the flat directions, and that the relic Q-ballSjWhich are stable against decay 
into both protons/neutrons (baryons) and neutrinos/electrons (leptons) |33[, may end up contributing to 
the energy density of dark matter [l^, [13]; thus, the Q-balls can provide the baryon-to-photon ratio [13], 
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i.e. n-hjn^ ~ (4.7 — 6.5) x 10"^" [s^ where n;, and are respectively the baryon and photon number 
densities in the Universe. 

In this paper we examine analytically and numerically the classical and absolute stability of Q-balls 
using flat potentials in the two specific models mentioned above. In order to study the possible existence 
of lower-dimensional Q-balls embedded in 3+1 dimensions, we will work in arbitrary spatial dimensions 
£); although of course the = 3 case is of more phenomenological interest. Previous work [2l|, [s^, [s^ 
on the gravity-mediated potential has used either a steplike or Gaussian ansatz to study the analytical 
properties of the thin and thick-wall Q-balls. Introducing more physically motivated ansatze, we will 
show that the thin-wall Q-balls can be quantum mechanically stable against decay into their own free 
particle quanta, that both thin and thick-wall Q-ball solutions obtained are classically stable against 
linear fluctuations, and conflrm that a Gaussian ansatz is a physically reasonable one for the thick- wall 
Q-ball. The one-dimensional Q-balls in the thin-wall limit are excluded from our analytical framework. 
The literature on Q-balls with gauge-mediated potentials has tended to use a test proflle in approximately 
flat potentials. We will present an exact profile for a generalised gau ge-mediated fiat potential, and show 
that we naturally recover results previously published in [H, [30l.l33| . 

The rest of this paper is organised as follows. In Sec. |TT] we briefiy review the important Q-ball 
properties that were established in [s^l . Section IIIII provides a detailed analyses for gravity-mediated 
potentials, and in Sec. IIVI we investigate the case of a generalised gauge-mediated potential. We confirm 
the validity of our analytical approximations with complete numerical Q-ball solutions in Sec. fVl before 
summarising in Sec. IVII Two appendices are included. In Appendix we obtain an exact solution 
for the case of a logarithmic potential, and in Appendix [Bl we confirm that the adoption of a Gaussian 
ansatz is appropriate for the thick- wall Q-ball found in the gravity-mediated potentials. 



II. THE BASICS 

Here, we review the basic properties of Q-balls as described in [s^l and introduce a powerful technique 
that enables us to find the charge Q and energy Eq of the Q-ball as well as the condition for its stability, 
and characteristic slope 7(w) = Eq/ujQ where to is defined through the Q-ball ansatz, which is given by 
decomposing a complex scalar field </> into = (T(r)e''^*. cr is a real scalar field, r is a radial coordinate, 
and therefore a; is a rotational frequency in the U(l) internal space. By scaling the radius r of the Q-ball 
ansatz, which minimises Eq, we can find the characteristic slopes in terms of the ratio between the surface 
energy S and the potential energy U of the Q-ball. When the characteristic slope, 7, is independent of 
w, we obtain the relation: Eq oc Q}^^ . In general the charge, energy and Euclidean action are given 

by 

Q = ^^ fr'; \^ (^'^" + ^-)' Eq=u,Q + S^, (1) 

where our metric is ds^ = —dt^ + hijdx^dx^ , the determinant of the spherically symmetric spatial 
metric hij is defined hy h = det{hij), and we have used the following notation: Jy^ = J d^x\/li — 

Hd-i /q°° dr r^^^, flo-i = v(d'/2) ' ^ ^^"-^ ^ number of spatial dimensions. Without loss of 

generality, we can take positive values of lo and Q. By defining the effective potential Ui^ of a potential 
U{a) 

U^ = U -^iv^a^ (2) 

the Q-ball equation is 

.^D-l, dU^ 

a H cr = -T—, 3) 

r d(T 

where a(r) is a monotonically decreasing function in terms of r. Given a potential U{a), which has a 
global minimum at cr = 0, it is possible to show that Q-balls exist within the restricted range of w Q: 

< uj < (4) 

where we have defined the lower limit tj^ = 2^1 , , > 0, cri (a;) is the nonzero value of cr where 

Ui^{(T+{uj)) is minimised (see Fig.[T|), and the upper limit lo\_ = . The existence condition Eq. ([4]) 

(T=0 
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restricts the allowed form of the potential 17, which implies that the potential should grow less quickly 
than the quadratic term (i.e. mass term) for small values of a. The case = corresponds to degenerate 
vacua potentials (DVPs), whilst uj- ^ has nondegenerate vacua (NDVPs). In [37| we examined the 
case of polynomial potentials and restricted ourselves to the case of uj\ = where m is a bare mass in 
the potentials. In this paper we extend our analysis allowing us to investigate the case ^ m^, needed 
since the potentials include one-loop corrections to the bare mass m. Here, the potential which wc will 
consider in the gravity- mediated models, is U = Ugrav + Unr where Unr is a nonrcnormalisablc term 
(to be discussed below), and 

Ugra. ^ \rr?cj-' {\ + K In ^ . (5) 

Here, if is a constant factor arising from the one-loop correction and i\/ is the rcnormalisation scale. 
To proceed with analytical arguments, we consider the two limiting values of lo ox = cr(r = 0) which 
describe 



I • thin- wall Q-balls when lo ~ a;_ or cquivalently uo ~ cr+('^)j ^g-j 
I • thick- wall Q-balls when uj ~ ujj^ or cquivalently (Tq — o'-(^)- 

Note, this limit doe not imply that a thick-wall Q-ball has to have a large thickness that is comparable to 
the size of the core size. For the extreme thin-wall limit, lo = ti>_, thin- wall Q-balls satisfy ~ ^{uj^)lo^. 
In particular, Coleman demonstrated that a steplike profile for Q-balls, which generally exist for uj- ^ 0, 
satisfies 7 = 1, which implies that the charge Q and energy Eq are proportional to the volume, and he 
called this Q-matter For absolutely stable Q-balls, the energy per unit charge is smaller than the rest 
mass TO for the field 4>, 

Thus the Q-ball satisfying Eq. ([7]) is stable against free-particle decays because the Q-ball energy Eq is 
less than a collection of Q free-particles of total energy Ej^ee = raQ. If the Q-ball has decay channels 
into other fundamental scalar particles that have the lowest mass nimin, we need to replace to by m^nin 
in the absolute stability condition Eq. ([7]). In the opposite limit uj ~ the Q-ball energy approaches 
the free particle energy, Eq — > toQ. For later convenience, we define two positive definite quantities, 
and by 



c l{c.'-u:l)al, (8) 
ml = n? - iJ- (9) 

which can be infinitesimally small for either thin- or thick- wall limits. By assuming {lo) ~ (i^-) = cr+ 
in the thin- wall limit, we immediately obtain the second line in Eq. ([8]). Notice that this assumption was 
implicitly imposed in our previous thin- wall analysis [37| . While this is fine for the gravity-mediated case, 
with Gauge mediated potentials which are extremely flat, this implicit assumption cannot hold because 
cr+(u;) does not exist. Therefore we will not use the variable e^^ for the case of the Gauge mediated 
potentials. Notice that the variable toJ cannot be infinitesimally small when we consider the gravity- 
mediated case: uP^ 7^ rr? . A powerful tool we can make use of when calculating some of the physical 
Q-ball parameters is the Legendre relation [13, [H] . For example the energy follows from 



Q 1^ ASu] 

>^UJ *■ V — ^ ~, 

diij 



Eq=ujQ + S^. (10) 



Assuming that 7 is not a function of uj, wc can compute the advertised characteristic slope, 

^=^^EQ^Q^h (11) 
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where we have used another Legendrc relation uj = 
classically stable, it satisfies 



in which we have fixed S^^. If a Q-ball is 



These classical stability conditions arc equivalent to the fission condition, i.e. ^ < in [s^ so that 
the charge Q for classically stable Q-balls is a decreasing function in terms of u>. By scaling a Q-ball 
solution with respect to the radius r, we also obtain the virial relation DU = —(£> — 2)S + DbjQ/2 and 
the characteristic slope 7(ti^), 

7H = 1+ (^-2 + i?|) (13) 

once the ratio S /U is given where S = Jy^ and U = Jy^ U are the surface and potential energies, 
respectively. For D > 2, we can see 7(w) > 1 because S, U > 0, which implies that is positive definite 
for D > 2, see Eq. ([T]), whilst S^i is positive for D = 1 only when U > S. It implies that we have to be 
careful to use the second relation of Eq. (fT2|) for D = 1 to evaluate the classical stability condition as 
we saw in the case of using the Gaussian ansatz, which is valid for D = 1 for polynomial potentials (37j . 
Our key results for Z3 > 2 are 

1 for 5 < U, 

(21?- 1)72(1?- 1) for5-W, (14) 
{D-l)/{D-2) for5>W. 

The first case in Eq. corresponds to the extreme case of thin and thick-wall Q-balls. Furthermore, 
in [s^l, we saw that for the extreme thin-wall Q-balls in DVPs, then there was a virialisation between 
S and which corresponds to the second case in Eq. ((M]). At present it is not known what kind of 
Q-ball potentials correspond to the third case; therefore, we will not be considering that case in the rest 
of our paper. Notice that in the case S for D = 2, we obtain the characteristic slope 7^1 from 
Eq. (fT5|) . Similarly, for D = 1, the characteristic slopes are obtained, i.e. 7—1, ^1, — 0, respectively 
for S ^ U, S ^ U, S ^ U. We will use these ID analytic results to interpret numerical results of 
one-dimcnsional Q-balls in the thin- wall limit. 

To end this section we note a nice duality that appears in Eqs. (|131 [H]) between the two cases S 
and S ^ Li. In particular, for S ^ U in D dimensions, the same result for 7 is obtained (to leading order) 
in 2 X Z) dimensions when S ^ U. 



III. GRAVITY-MEDIATED POTENTIALS 

The MSSM consists of a number of flat directions where SUSY is not broken. Those flat directions 
are, however, lifted by gauge, gravity, and/or nonrenormalisable interactions. In what follows the gravity 
interaction is included perturbatively via the one-loop corrections for the bare mass m in Eq. ([5]) and 
the nonrenormalisable interactions (Unh), which are suppressed by high energy scales such as the grand 
unified theory scale Mjj ~ 10^^ GeV or Planck scale iripi ~ 10^^ GeV. Here, m is of order of SUSY 
breaking scale which could be the gravitino mass ^ 7773/2, evaluated at the rcnormalisation scale M (23| . 
We note that, following the majority of work in this field, we will ignore A-term contributions ( U(l) 
violation terms), thermal effects [s^ H^l which come from the interactions between the AD field and 
the decay products of the inflaton, and the Hubblc-induced terms which gives a negative mass-squared 
contribution during inflation. It is possible that their inclusion could well change the results of the 
following analysis. 

The scalar potential we arc considering at present is [U, [2^ 

U = Ugrav + UnR = ^777^^2 ( 1 + K In ( ^)) + ^^'t'' (15) 



\ \ / / pi 

where we used Eq. ([5]), if is a factor for the gaugino correction, which depends on the flat directions, and 
M is the rcnormalisation scale. Also A is a dimensionless coupling constant, and Unr = n-i c", where 
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n > 2. If the MSSM flat directions include a large top quark, K can be positive and then Q-balls do not 
exist. For flat directions that do not have a large top quark component, we typical find K ~ —[0.01 — 0.1] 
[21I . I41I I . The power n of the nonrenormalisable term depends on the flat directions we are choosing along 
which we maintain R parity. As examples of the directions involving squarks, the u'^dfiF direction has 
n = 6, whilst the u'^u'^<Fef direction requires n = 10. A complete list of the MSSM flat directions can 
be found in Table 1 of [4^. Since the potential in Eq. (fT5|) for K <0 could satisfy the Q-ball existence 
condition in Eq. (j4]), where w+ 3> m, Q-balls naturally exist. 

In the rest of this paper, we will focus on potentials of the form of Eq. ([TS]) for general D{> 1) and 
Lo and n{> 2) so that M and m^; have the same mass dimension, {D — l)/2, as a. It means that the 
parameters M and rupi are only physical for D — i. For several cases of n and D, the term C/^rfl can be 
renormalisable, but we will generally call it the nonrenormalisable term for the future convenience. The 
readers should note that the potential Eq. has been derived only with M =1 supcrgravity in D = 3; 
therefore, the potential form could well be changed in other dimensions. Furthermore, the logarithmic 
correction breaks down for small a and the curvature of Eq. p5|) at c = is flnite due to the gaugino 
mass, which affects our thick- wall analysis and their dynamics. However, we concentrate our analysis 
on this potential form for arbitrary D, n and any values of a for two main reasons. The flrst is that it 
contains a number of general semiclassical features expected of all the potentials, and the second is that 
it offers the opportunity to consider the lower-dimensional Q-balls embedded in D = 3. 

In Appendix 13 wc obtain the exact solution of Eq. ^ with the potential U = Ugrav', however, exact 
solutions of the general potential U in Eq. psp arc fully nonlinear and can be obtained only numerically. 
Therefore, wc will analytically examine the approximate solutions in both the thin and thick-wall limits. 
Before doing so, we shall begin by imposing a restriction on A in Eq. psp in order to obtain stable 
Q-matter in NDVPs. With the further restrictions on A and wc can proceed with our analytical 
arguments, and we will finally obtain the asymptotic Q-ball profile for large r which will be used in the 
numerical section, Sec. [V] 



A. The existence of absolutely stable Q-matter 



As we have seen, the first restriction on the gravity-mediated potential Eq. (|15p which will allow for 
the existence of a Q-ball solution Eq. ^ is K < 0. However, given values for m, nipi, M, n, and K in 
Eq. psp . we need to restrict the allowed values of the parameter, A in the potential in order to ensure we 
obtain absolutely stable Q-matter. Notice that Q-matter exists in NDVPs, whilst the extreme thin-wall 
Q-balls in DVPs, which will not be Q-matter as it will turn out, may exist with the lowest possible limit 
of A. 



By using the definitions of u)^ and (T+ , namely, uj'^ 



m and 



da 



0, we shall find the range 



of values of A for which absolutely stable Q-matter solutions exist. Moreover, we will obtain the curvature 
/i, which is proportional to \K\, of the effective potential at cr+. 

The effective potential for Eq. (|15p can be rewritten in terms of new dimensionless variables = 

cr /M, u) = Lj/m, and 



|ApM"-^ 

Tl — 4 9 



>o, 



(16) 



1 



= ^M'^m'a^ (1 - - 2\K\ Ina) + A/^m^/J^a" 



(17) 



After some simple algebra and introducing Co^ = ^ 
li^l 



{n - 2)/32 



and 



0, we obtain 



1 



n - 2 



71 -2 + 21X1 -21X1 In 



\K\ 



{n - 2)/32 



(18) 



Notice that Cj'^ = corresponds to DVPs where Q-matter solutions do not exist [Sjl, whilst the extreme 
thin-wall Q-balls do exist and are absolutely stable as wc will see. In NDVPs, Q-matter solutions exist 
and arc absolutely stable when < Cu'^ < 1, see Eq. ([7]). Combining these facts and using the second 
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relation in Eq. (jlSp . we have the constraint on A for stable Q-matter solutions to exist, namely 



\K\e- 



\K\e 



n 

n — 



■ cxp 



n ■ 



-1 to", ^m? 



2 M 



n-2 



■ exp 



2\K\ 
n-2 



< < 



< lAI- 



< 



\K\e- 



n-2 M 



n-2 



(19) 
(20) 



where we have used Eq. to go from Eq. to Eq. Here, the lower limit of |Ap corresponds 

to o)^ 0, whilst the upper limit corresponds to oj'^ ~ 1. The inequality in Eq. (|20p implies that if 
the coupling constant A of the nonrenormalisable term in Eq. (jl5p is too small, then it docs not support 
the existence of Q-balls. whereas a large A coupling leads to unstable Q-mattcr. With the following 
parameter set, m = M = 1, \K\ — 0.1, n 6 and the lower/upper hmits of in Eq. p^ . Fig. [T]shows 
the inverse potentials in Eq. p7p and their inverse effective potentials —U^^ with various values of w. The 

lower limit, 0^ — ^^^^ 



0" = 



■exp (^^i^j, corresponds to DVPs case with = 0, whilst in the upper limit, 

^ , the potentials do not have degenerate vacua with uj^ — 1, hence are called NDVPs. By 
substituting the values of into Eq. (US]), we obtain the values of (t+ indicated in Fig. [TJ Finally we can 

,21. ,\ — d^C/^ 



obtain the curvature, {to) 



, evaluated at ti>_ , i.e. 
2 = /i2(cj_) = m^\K\{n - 2) oc 



(21) 



which implies that a small logarithmic correction \K\ ^ in Eq. gives an "extremely" flat 

effective potential compared to the quadratic term w? around a = a+ for a given n ^ 0(10°^^). 
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FIG. 1: Parameters a±{Lu) for a potential of tlie form U{a) — ^a^ (l — \K\ Ina^) +/3^(7^ (effective potential Uu = 
U— ^uj^a'^) with li^l = 0.1. The left hand figure corresponds to the case of a DVP with = ^-^^ — exp ^ — ~ 

1.90 X 10"", whilst the right hand side is the NDVP with (3'^ = ^^^f^ ~ 9.20 x 10"^, see Eq. pH]). The coloured 
lines in each plot correspond to different values of uj. The variable ct+{u) is defined as the maximum of the inverse 
effective potential —Uuj whereas a-{Lj) corresponds to — !7ij(o--(tj)) — for cr_(aj) 7^ 0. Recalling uj- = in DVP, 

the DVP has degenerate vacua at o-+{0) = e^''^ exp ^j^^ ~ 1-91 x 10^ (red-solid line), whilst the NDVP does 
not. The inverse effective potential —Uui with uj- = 1 in NDVP (green-dashed line), however, has degenerate 
vacua at cr+(cj-) = e^^* ~ 1.28, see the first relation in Eq. (|18|l . For the lower limit lj ~ uj- (green-dashed 
lines), we could see a+ — e^^^, whilst the purple dotted-dashed lines show (y-{uj) near the thick-wall like 
limit Lo = 3.0 ~ ui+ where ui+ S> 1. 



B. Thin-wall Q-ball for ao ~ a+, Rq » S, l//i, D>2 

For the extreme limit ui = uj-, Coleman demonstrated that the steplike ansatz 0] is applicable to 
the case of NDVPs because the surface effects of the thin- wall Q-ball in this limit are not significant. 
There arc situations though where we would like to explore the region around lu = uj^, corresponding 
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to ao — (J+(cj), and to do this we need to include surface effects. In [33| we explained how to do this 
under the assumptions: Rq/S, fiRq ^ 1, ct^Rq) < cr_(u;), a+{ui) ~ (T+(a;_) = (t+, and that the surface 
tension r ~ * du^JHTZZ does not depend "sensitively" on w. Here, Rq, 5 are, respectively, the Q-ball 
core size and the shell thickness. We note that in [4^, Coleman assumed C/^ ~ Ui^_ in the shell region, 
and this is equivalent to saying <j{Rq) < cr-(w). In what follows we will be making use of Coleman's 
approach. By requiring this or (j{Rq) < a-{uj), we can guarantee real values of shell thickness 5 and 
surface tension t. The assumption, in which r does not depend on w, is related to the assumptions: 
cr+{uj) ~ (T+ and 11^; — Uuj_ is negligible in the shell region. 

Under these assumptions and for Z? > 2, we now apply the previous thin- wall analysis in [37| to the 
present potential Eq. (jlSp . The ansatz is given by 

{cr+ - s{r) for < r < Rq, 
a{r) for RQ<r <Rq+5, (22) 

ioT Rq+5 < r, 

where Rq, S, the core profile s{r), and the shell profile a{r) will be obtained in terms of the underlying 
potential by extremising Su, with respect to Rq. Each of the profile hmctions satisfies 

" ^^^-s' - ^is = 0, (23) 



r 
a" 



dU.., 



da 



(24) 



By recalling Eq. ([8]), we have previously found that [37 



Rq ~ {D-l)—; S^~-dVD>0; Q~ujaiVD, (25) 
Eq 2D - I ujI 

ZQ ~ 2{D-1) ^ 2(Z?-l)a;2' ^^^^ 

oj dQ 2L>w2 ^ ^ 

nTT ^ ^<0' (27) 

Q duj UJ — 

where we have taken the thin-wall limit lo ~ w_ in the last inequality. Notice that our analytical work 
cannot apply for the \D thin- wall Q-ball, sec the first expression in Eq. (|25p . 

NDVPs: This type of potential supports the existence of Q-matter that corresponds to the regime 
U ^ S. The Q-matter can be absolutely as well as classically stable for the extreme limit w ~ when 
the coupling constant A for the nonrcnormalisable term in Eq. (|15p satisfies Eq. ((20|) . The characteristic 
slope is given by the first case of Eq. p4p . and the charge and energy are linearly proportional to the 
volume Vd- 

DVPs: With the presence of degenerate minima in Eq. (fTSj) . in [s^l we obtained the ratio U/S ^ 1, 
which corresponds to the second case of Eq. (fT4| . The charge and energy are not proportional to the 
volume Vd itself in this case; hence, we cannot sec the existence of Q-matter in the extreme limit 
w = w_ = 0. Instead we can find the proportional relation simply from Eq. and Eq. (|26p . namely 

SqCxQ2(C-1)/(2C-1), 

Our main approximations are based on the assumptions ctq ~ cr_|_, Rq ^6, and ~ Uuj_ in the 
shell region. In what follows we will see through numerical simulations that our analytic results agree well 
with the corresponding numerical results even in a "flat" potential choice \K\ = 0.1, m = M = 1, n = 6, 
which implies that 1/fJ. ^ 1.58, see Eq. (PT|) . 



C. Thick-wall Q-balls for /^^ < \K\ < 0(1) 

In [s?! we studied thick- wall Q-balls in general polynomial potentials, and extracted out the explicit 
u> dependence from the integral in S^^ by reparameterising terms in the Euclidean action S^, in terms of 
dimensionless quantities and by neglecting higher order terms. We then made use of the technique Eq. PH)) 
and obtained consistent classical and absolute stability conditions, Eqs. ([71 [T2)) . For our present potential, 
Eq. ([T7)) . which satisfies the condition, < \K\ < 0(1), wc will be able to ignore the nonrcnormalisable 
term by introducing a = a/M and in Eq. (jl6p . Wc can then obtain the stability conditions using the 
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same technique as before. Indeed for the Hmit oj > 0{m), we will see d{r) ~ 0{e) < 0(1) where e is a 
small dimensionless constant (not e^} in Eq- (IHl))j and see ctq = i5'(0) > d{r) for any r because d{r) is a 
monotically decreasing function in terms of r. Since the leading order of the logarithmic term, CT^lna, 
in Eq. p7|) is of O(e^) using the L'Hopital's rules, we can ignore the nonrenormalisable term in Eq. (fT7|) 
at the beginning of our analysis. To confirm this, in Appendix [B] wc will keep all terms in Eq. (jl7p by 
introducing a Gaussian ansatz and show that the results below [Eqs. ((33l [34]) ] can also be recovered under 
the same assumption < \K\ < 0(1). By adapting the techniques introduced in [s^], in this subsection 
we will show how to obtain the thick- wall solutions without involving the Gaussian ansatz. 

First of all we introduce two characteristic limits: the "moderate limit" lo > 0(m) and the "extreme" 
limit w ^ m. We will see (Jq ~ d'-{uj) O"*" which leads to a'_(<^) ^ 0(1) in the "extreme limit", and 
then even in the "moderate limit" we will sec that the contributions from the nonrenormalisable term 
are negligible and that (T_(a;) is a monotonically decreasing function in terms of uj. Under the conditions 
^ \K\ < 0(1) in Eq. we obtain 



(28) 
(29) 
(30) 



\K\w? da-{uj) 



-1 + 2{n - 2) 



\K\ 



2\K\tv? 



-1 < 0, 



0, 



where we used Uu){cr-{uj)) = to obtain Eq. ([28]). It follows that cr-{uj) ^ 0(1) for the thick-wall 
limit a; ^ TO, and wc can ignore the nonrenormalisable term. Since Eq. (j29p implies that '^"'^J"-* < 
in the limit o^{uj) < 0(1), ^-{io) is a monotonically decreasing function. Therefore, we can ignore the 
contributions from the nonrenormalisable term up to w > 0(to) which wc call the "moderate limit" with 
the notion as seen in the second relations of Eqs. ((28l [29]) . instead of the "extreme" limit S> m with 
the notion Thus, we obtain the desired results of the second relation in Eq. ([30]) . From Eq. ((28|) . the 
logarithmic term may be of < 0(1) for \K\ < 0(1), ^ 0(1) in the "moderate" limit, which implies 
that the "moderate limit" is valid even when lo ^ 0(to). 

Let us define a{r) and r through (T(r) = aa{r) and r = br where a and b will be obtained in terms 
of the underlying parameters. By substituting these reparamerised parameters a, f and neglecting the 
nonrenormalisable term into Eq. ([T]) due to 'the L'Hopital's rules', we obtain 




- 2\K\lna 



= a^M^b 



2iD-2 



S{a), 



(31) 
(32) 



where S{a{f/b)) = ilo-i J drr^ ^ |i (^)^ — 5^1^(1 — lna^)|, which is independent of lu. In going 

from Eq. pT|) to Eq. ([5^ we have set the coefficients of the three terms in the brackets of Eq. (PT|) to be 
unity in order to explicitly remove the lu dependence from the integral in S^. In other words, we have 



set a 



-1/2 



exp 



2\K\7n^ 



'^/^a^iuj), b - 



Following Eq. ([TU]), we can differentiate Eq. ([32 



with respect to lu to obtain Q and then use the Legendre transformation to obtain Eq. Coupled with 
Eqs. ((7l fT2)) we obtain both the classical and absolute stability conditions. This is straightforward and 
yields 



2lu 



_d_ fEQ 
duj \ Q 



m^\K\ "' 



1 - 



2lu^ 



Eq m^\K\ 
tuQ^ 2^2 

-.1>0, ^^-1- 

Q dhj 



2uj^ 



2u? 



rn^lKl 



<0, 



(33) 
(34) 



where we have taken the "extreme" limit lu ^ m as indicated by '— >'. Equation (j33p implies that the 
characteristic slope for the thick- wall Q-balls are tending towards the case S -^U in Eq. and Eq. 
shows that the Q-balls are classically stable. These results are independent of D. In Appendix IB] we 



9 



will generalise the results of Eqs. ([551 [M]) by adopting an explicit Gaussian ansatz without neglecting the 
nonrcnormalisable term. 

Before finishing this subsection, let us comment on possibilities to have absolutely stable thick-wall 
Q-balls in the case, \K\ < 0(1), 0^ ^ 0(1). The results present above still hold even in the "moderate 
limit" Lo ~ 0{m) for the present case. Thus, the thick-wall Q-balls, if they exist, can be absolutely stable 
when the following conditions from Eqs. (|7l I33p are met: 



UJ 1 + ^l-2\K\ , ^, 1 ^2 



^-<m, -< ^- ^, < -, « 0(1). (35) 

TO 2 2 

It follows that for \K\ > 1/2, the thick-wall Q-balls are always absolutely unstable. If ui- > to, we know 
w > UJ- in both the "moderate" and "extreme" limits, hence the thick-wall Q-ball is always absolutely 
unstable again, see Eq. ([7]). Notice that the condition < \K\ implies uj- < 0{m), see Eq. p^ . so 
the first condition in Eq. (|35[) can be satisfied. This then leaves only a small window of the parameter 
space for absolutely stable thick-wall Q-balls. In the numerical section. Sec. |Vl we will confirm that 
the thick-wall Q-ball can be absolutely stable against decay into their own quanta by choosing suitable 
parameters, i.e. uj- =0, ~ 1.90 x 10"^-^, and \K\ = 0.1. 



D. Asymptotic profile for large r and (3^ < \K\ < 0(1) 

In order to obtain the full numerical profiles over all values of uj, wc should analytically determine the 
asymptotic profile for large r in the potential Eq. ^ which satisfies <\K\<0{1) as in the previous 
subsection. As long as the value of r satisfies r > where R^u is some large length scale and depends 
on UJ, we can assume that the friction term in Eq. Q and the nonrcnormalisable term in Eq. (|15|) are 
negligible for large r. Hence, the Q-ball equation Eq. ([3]) reduces to the one-dimensional and integrable 
form 



da 



(36) 



where ~ ^rn^a^ (^1 — ^ \-^\ ^'-'S ("^j j ' -'^l^^^i'^^ ([55]) implies that the profile has a symmetry 

under the variation of r because Eq. (j36p does not depend on r explicitly. Multiplying both sides of 
Eq. (EH) by ^ leads to 

"'^'^'^ da 

= R^-r, (37) 

where we have used the boundary conditions: a'{oo) 0, Uuj{a{oo) 0) ^ and a'{r) < 0. After 
some elementary algebra, the final asymptotic profile becomes 



a{r) = Me^^</2™^ exp ( _!!!!^(, _ ,^)^ ) , (33) 



d_ (a^ 
dr \ a 



= m^\K\M, (39) 



where r^^ = R^ — ^ ^ — '^^^ log (^^^^7^^ /{\K\m). Equation ([55]) is a consequence of the symmetry in 

Eq. (j36p under the translation r ^ r ~ from a Gaussian profile as seen in Eq. (jAip of Appendix lAl 
Furthermore, Eq. ([39)l depends on the parameters to, M, \K\ in Eq. We will later use the relation 

Eq. (|39p as a criterion that must be satisfied in obtaining full numerical profiles for all values of uj. 

Wc finish this section by recapping the key results we have derived for the case of the gravity-mediated 
potential, Eq. (jl5p . in both the thin and thick- wall limits. In the thick-wall limit, wc imposed the 
restrictions 0^ < \K\ < 0(1) on the potential to ignore the nonrcnormalisable term. In both limits, we 
have derived the characteristic slopes in Eqs. (PSI [55)1 and the classical stability conditions in Eqs. (^71 IM)) 
and shown that the Q balls are classically stable in both cases. The thin- wall Q-balls in DVPs are always 
absolutely stable, and Q-matter in NDVPs can be absolutely stable when the coupling constant for the 
nonrcnormalisable term satisfies Eq. (|20p : whilst absolutely stable Q-balls in the thick- wall limit may 
exist only for Eq. ([551) . Finally, we obtained the general asymptotic profile, Eq. ([SS]) . for large r. 
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IV. GAUGE-MEDIATED POTENTIAL 



The Gauge mediated scalar potential can be written in quadratic form in the low energy regime for 
scales up to the messenger scale Ms, and carries a logarithmically (extremely) flat piece in the high 
energy regime [l], [l^] ■ This extreme flatness means that the thin- wall Q-ball we used in Eq. cannot 
be applied to this situation, and so we now turn our attention to Q-balls in extreme flat potentials. We 
will generalise the results of [l|] to an arbitrary number of spatial dimensions and show that the known 
Q-ball profiles in P.i33| are naturally recovered by our more general ansatz. Moreover we will investigate 
both the classical and absolute stability of these Q-balls. The gauge-mediated potential, which we will 
use in this section, is approximated by [13, 

U{a)^[}f^" ^ M'Jr/^/^i' (40) 
^ ' — const, for cr(i?) < air), ^ ' 

where C/q and R are free parameters that will be determined by imposing a condition that leads to a 
smooth matching of the profiles at cr(i?), Uq ~ ^m^a^{R). Notice that Q-balls exist within < uj < m 
in Eq. (HO]), and the potential does not have degenerate vacua although uj^ ~ 0. Since Eq. is not 
differentiable at (j{R), we can approximate Eq. (PD|) by 

U„ = lm^A'{l-e-^'/^') (41) 

which we will use in the numerical section. Sec. |Vl Note that A = cr(i?) corresponds to the scale below 
which SUSY is broken, so that Uq = ^m^A^ in Eq. ([40]). The potential Eq. ([4T]) differs from the one 
used in [i^ , but is similar to the potential used in [43] . Fig. [2] shows the inverse potential Eq. (|^T|) 
and the inverse effective potentials for various values of w with m = 1, A^ = 2, which implies Uq = 1. 
The red-solid line shows the inverse potential of Eq. ((4T|) {—U gauge) , and the sky-blue dotted-dashed 
line corresponds to the inverse quadratic potential of Eq. (j40|) . For sufficiently large and small cr, the 
two potentials in Eqs. (j40i I4ip have similar behaviour, but we can see the difference in the intermediate 
region of a where 1 ^ cr < 3. Hence, we can expect that profiles around the thick- wall limit are different 
between the potentials since the thick-wall profiles are constructed in the particular region, 1 ^ cr < 3; 
hence it may lead to the different stationary properties and stability conditions. 

Using Eq. ([8]), the Q-ball equation, Eq. ([3]), in the linearised potential Eq. ((40|) becomes 

encore + — ^O'core + ^'^(^core = 0, for < r < i?, (42) 

r 

<^sheii + — -<^'sheii - mlasheii = 0, for i? < r, (43) 

r 

where the profiles should be imposed to satisfy the boundary conditions, a' < 0, cr(0) = (Tq = 
finite, a{oo) — a'{oo) = 0, cr'(O) = 0. The solutions are 



a,ore{r) = A r^-°I^JDIi-i{^r) for < r < i?, 
Osheiiir) = B for R < r, 



(44) 



where J and K are Bessel and modified Bessel functions respectively, with constants A and B. By 
introducing ctq: and expanding Jd/2-i(w) for small cor in acoreir), and by using the condition Uq = 
I'^^'^^sheiii^) we obtain 



UJ 



D/2-1 



1 



A = aoT{D/2)(-\ , Uq = -m'B'R'-'' Kj,,^_^{m^R). (45) 



Since the energy density is smooth and finite everywhere, we have to impose a smooth continuity condition 
to the profiles crcore{R) = cTsheiiiR) and u'coreiB) = Kheiii^)^ which from Eq. dH]) gives 

A _ Ko/2~i{nT'ujR) _ m^KD/2{m^R) 



We will see that the particular value of ctq does not change important features such as the stability con- 
dition and characteristic slope of the Q-ball solutions. Using Eq. (|46|) we obtain the following important 
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FIG. 2: The inverse potential —Ugauge in Eq. (|4ip (red-solid line) with m = 1, = 2 which implies f/o = 1 
and the inverse effective potentials —Uuj for different values of uj. In order to compare between Eq. (|40p and 
Eq. (|4H) . we plot the inverse quadratic potential with the sky-blue dotted-dashed line. The two potentials are 
asymptotically similar, but they are different around the intermediate region of a, where 1 < o- < 3. 



identities, which we will make use of later [45 



LU- 



Jd/2-i{^R) 
Jd/2{^R)Jd/2-2{^R) 



KD/2{m^R) 
KD/2-i{mujR) ' 

m^\2 Ku/2{mu:R)KD/2-2{mu,R) 
LU / 



^D/2 



where we used the recursion relations Jp_i(z) + ,7^-1-1(2:) = ^Jf^{z), — = —^Kf^{z) 



for any real 11 and z. We can easily find 
- {\(7sheii{R)(y'sheiiiR)) '^^D , and then using 

5, 



UoVd 



^a,ore{RycoreiR)) dVo and 5^'-" 



(47) 
(48) 



gcore ^ gshell foWo^s that 

(7ol^D, (49) 

where we have again used the continuity relations (Jcore{R) = crsheii{R) and CT^Q^g(i?) = ^'sheiiiR)- 
find the charge Q, we do not make use of the Legendre relation Q — — in Eq. pO|) . because i? is a 
function of w, and is determined by Eq. ([17|) . However, we can obtain Q by substituting Eq. directly 
into Eq. Q: 



DUoVd ( A'c/2(m^i?)A'£,/2-2(w<^^) 



^D/2- 



(50) 



where we have used Eqs. (^51 and Eq. ([35]), as well as the relation, 

/ dy yZf^{y) = ^ (^^(y) — Z^-i{y)Z^^i{y)) , [H, |4l|. Here, /i is real, and Z can be either the Bessel 
function J or the modified Bessel function K, and we have used the following recursion relations to obtain 



±z J, 



the indefinite integral: 2: -37- ± /iJp = =nz J^=pi, J^-i 

For future reference we obtain explicit expressions for R for case with an odd number of spatial 
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dimensions. Eq. ([17|) can be solved explicitly in terms of R to give 



ujR = arctan ) , for Z) = 1, (51) 



ujR = TT — arctan ) , for D = 3, (52) 



where we have used J3/2{x) = (^^T^ ~ cos(a;)^ , Ji/2{x) = y^sin(.T), J-i/2{x) = 

cos(a;), K3/2{x) = ^/^e^'' (l + i) , Ki/2{x) = ■f^-i/2(a;)- We will discuss the classical 

stability for Q-balls in D = 1, 3 in the numerical section, in which we will show stability plots arising 
from Eqs. dSmnH). 

A. Thin-wall-like limit for m^R,ujR > 0(1) 

We now discuss both the classical and absolute stability of gauge-mediated Q-balls in arbitrary dimen- 
sions D, in the limit nit^R, ujR ^ 1, which implies that the "core" size R is large compared to l/m^^, l/u. 
As we will sec in the numerical section, Sec. |Vl the limit will turn out to be equivalent to the thin-wall 
limit LU ~ ~ 0. Recall that this potential docs not have degenerate vacua. Using Eqs. P51 ISO)) . 

S^^^{l + 0{im^Rr')}, (53) 
where we have used linii^i^oo ^^{^) ^ ^ + + 0{z~'^) 



Ea D + l 



The characteristic slope follows 

.Q- D (^4) 

from which we see immediately from Eq. pT|) that we recover the published results of [l|, , namely 
E (X From Eqs. QEH), the thin-wall-like Q-ball is absolutely stable since the present limits 

will cover the thin-wall limit uj ~ lj- ~ as we stated. 

We can also obtain an explicit expression for R{uj) and ^ in the limits m^R 3> 1 and loR ^ ^ 1:1, 
where /i is the argimrcnt of the Bessel function: 

uR = ( — i — ) TT — arctan [ ), (55) 



4 / \m, 

= 1 ■ (56, 

dio ^ \ m^R / Lo 

Notice that Eq. for D = 3 reproduces the given profile in P, [l^ , and it coincides with the exact 
expression derived in Eq. ((52|) . Using Eqs. ((50l [55]) and Eq. ((56)) . we obtain 



VdUoD 



(57) 



^^--i^^-KO, (58) 

which shows that the Q-ball in this limit is classically stable. One can also check both Q ~ ^lif = 
DUqVd/uj from Eq. ((56)) and ^ ("^) — ^^W' ^ from Eq. ([54]) . which are respectively consistent with 
Eq. dnSl) and with the resuh in Eq. (fH)) . 

B. Thick- wall limit for D = 1, 3, ... 

Having just discussed the thin-wall-like properties for arbitrary D, we turn our attention now to the 
the other limit, uj ~ uj^. This is much more difficult to analytically explore because Eq. (|48[) can only give 
a closed form expression for R for the case where D is an odd number of spatial dimensions. Therefore, 
we will concentrate here on the cases, e.g. D = 1, 3. 
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D = 3 case: From Eq. (|5^ and recalling that in the thick-wall limit, to^^ ^ 0, w ~ uj+ = m, we 

dF 



obtain i? ~ ^ ^ g^^j ]jy substituting these into Eq. (|50|) we find 



uj dQ I 
Q duj 



-1 + — ^ — >0, (59) 



^ = l + -r^-l, (60) 

which shows that the three-dimensional thick-wall Q-ball is classically unstable. This fact is consistent 
with the relation that (^) = 1 - < where we have used Eq. It also follows 

that the thick-wall Q-ball is not absolutely stable, and the solution will decay to free particles satisfying 
Eq mQ which is the first case of Eq. 

D = 1 case: As in the case Z) = 3, Eq. ([5l|) implies i? ^ 0, ^ ^ ~;;r^ thick- wall limit. 

Using the above results, we obtain 



UJ dQ ^ rn? 

Q duj ijp' 

J^Q _ , (- 1 



>0, (61) 
L. (62) 



Note that the approximate value in Eq. (|6ip is the same as Eq. ((59| . Then the one-dimensional thick- 
wall Q-ball is also classically unstable. This fact is again consistent with the result that ^ (■^) — 

1 + rUujR ~ ^ — —iiLJi < 0. As in the three-dimensional case, the thick- wall Q-ball is not 

absolutely stable, and the solution decays into its free particles. 



C. Asymptotic profile 



The asymptotic profile for the large r regime in this model can be described by the contribution from 
the quadratic term in the potential Eq. (j40p . from which the profile is 

<r) ^ E.[^r-^e-"'-'- ^ -- ^ + (63) 

V ZttIui (J 2r 

where i? is a constant (33] • Note that we have used the fact that the modified Bessel function of the 
second kind has the relation K^{r) ~ \f^^~'^ large r and any real /x. We will use the criterion in the 
second expression of Eq. ([55)1 in the following section. 



Summarising our most important results, the thin- wall-like Q-ball is classically stable for a general 
Z), whilst it is absolutely stable as seen in Eqs. (|541 [58]) . On the other hand, for thick- wall Q-balls in 
D = 1, 3, the Q-balls are both classically and absolutely unstable, as can be seen from Eqs. (|59l [60|) and 
Eqs. ([61 1 [62 ]) . Finally we obtained the general asymptotic profile Eq. ([63]) for large r. 



V. NUMERICAL RESULTS 



In this section, we obtain exact numerical solutions for Q-balls for both the gravity-mediated potential 
in Eq. (fT7|l and the gauge-mediated potential in Eq. (|4T]l with dimensionless parameters by setting m = 
M — \ and A^ = 2. We adopt the 4th-order Runge-Kutta algorithm and usual shooting methods 
to solve the second order differential equations Eq. ([3]) (for full details see the numerical techniques 
developed in (37|V The raw numerical data contains errors for large r, thus we introduce the previously 
obtained analytical asymptotic profiles to help control these uncertainties. In particular we use Eq. p9p 
for the gravity-mediated potential and Eq. for the gauge-mediated case. Using these techniques, 
the numerical profiles match smoothly and continuously onto the analytic ones. In order to check the 
previously obtained analytic results, we calculate Q-ball properties numerically over the whole parameter 
space Lo except around the extreme thin- wall limit w = w_, because it is difficult to obtain reliable 
numerical results in that limit. 
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A. Gravity-mediated potential 

Wc shall investigate the gravity- mediated potential with two choices of A in Eq. ^TE\i for \K\ = 0.1 and 
n = 6, which can be seen as the red solid lines in Fig.[TJ The choice of the parameters, \K\ and n, are 
simply from phenomenological reasons. The degenerate vacua potential (DVP) on the left has (jJ_ = 

(/32 = i^^^exp (-^) ~ 1-90 X 10"" < and the nondegenerate vacua potential (NDVP) on 

the right has cj_ = 1 (/S^ = IBfl ^ 9.20 x 10"^ ^ 0(1)), recalling Eq. (HH). Fig. [1] also shows plots of 
the inverse effective potentials —U^ for various values of tu. Because of numerical complications, we are 
unable to fully examine the properties in the extreme thin-wall limit; however, by solving close to this 
wall limit, our numerical results recover the expected analytical results we derived in Eqs. ([26l [27|) . With 
the above choice of parameters, the curvature fi of at (J-^{uj-) = a+ in Eq. (|2ip is /i^ ~ 0.4 which 
implies that l//i ^ 1.58. From the first relation in Eq. pS)) . we have found a^^ ^ 1.28 in NDVP and 
(T+ ^ 1.91 X 10^ in DVP. Since we have assumed Rq ^ l/fi, cto — 0-+ in our thin- wall analysis for the 
gravity-mediated potential, we see that it breaks down when the core size Rq becomes the same order 
as 1/^ and/or ctq 7^ cr+. Although the full definition of the core size Rq is presented in [sH, it is very 
time consuming to evaluate it properly in the simulations; hence, in this analysis we have used a more 
naive approach, in which we have estimated the value of r = Rq when the field profile drops quickly from 
its core value. For the thick-wall limit, we required the condition < \K\ < 0(1), which is satisfied 
with the above chosen parameter set; hence, the analysis is valid for oj > 0(1). Because of the choice of 
\K\ = 0.1 < 0(1) and uj- = in NDVP, we will see our analysis holds even for oj ^ 0(1). 

Hybrid profile: The numerical profiles have errors for large r which correspond to either undershoot- 
ing or overshooting cases; thus, to minimise the errors in the region of large r we replace the numerical 
data by the predicted asymptotic analytical profile using the criterion Eq. (|39p to obtain the solution for 
the whole range of r. We then have the hybrid profile which can be written as 

for r <C Rnum : 

(64) 

for /? < r < R 

where (Tnuin is the numerical raw data, Rnum is determined by | {^o'^^^j Onunxi — l|r=_R„„„, < 0.001, and 
we have set Rmax = 60 throughout our numerical simulations in this subsection. We have calculated the 
following numerical properties using the above hybrid profile, Eq. (|64p . for D = 1, 2, 3: 

Profile: In the top two panels of Fig. [3] (DVP on the left and NDVP on the right), the red-solid and 
blue-dotted lines show the numerical slopes —a'ja for two typical values of w in I? = 3. We smoothly 
continue them to the corresponding analytic profiles by the methods just described in the numerical 
techniques, see green-dashed and purple-dotted-dashed lines. The linear lines correspond to the Gaussian 
tails in Eq. and for the cases of w = 0.14 (DVP) and cj — 1.01 (NDVP) corresponding to the thin-wall 
solution we see that it is shifted from the origin to r ~ 21. The middle panels show the obtained hybrid 
profiles of Eq. ([M)) for the various values of uj and D. The higher the spatial dimension, the larger the 
core size Q-balls can have. The energy density configurations pE^f) can be seen in the bottom panels of 
Fig. m Outside of the cores of the DVP profiles for cj ~ uj- , we can see the same features as we saw in 
the polynomial potentials we investigated in (37l |. namely, highly concentrated energy density spikes. In 
NDVP, however, the spikes cannot be seen. The presence of the spike contributes to the increase in the 
surface energy 5, which in turn leads to the different virialisation ratio for S jU where U is the potential 
energy, as can be seen in Eq. p^ . 

Criterion for the existence of a thin-wall Q-ball: Fig. |4] shows the numerical results for aQ(yj) 
against u) for both types of potentials - DVP (left) and NDVP (right). Our main analytical approximation 
relies on gq(uj) ~ (y+{yj) ~ (t+ = aj^[u)-\ where we have found Oj^ ^ 1.28 ~ 0(1) in NDVP and 
cr+ ~ 1.91 X 10^ S> 0(1) in DVP. The 3D thin-wall Q-ball (green-crossed dots) appears for a wider range 
of u) than the ID Q-ball (red-plus signs) in DVP as well as NDVP. For each case, the approximation 
can be valid, respectively, up to w '--^ 0.24 or w ~ 1.04 with about 10% errors for the ZD case. Near 
the thick- wall limit uj ~ ujj^ for both potentials, we see (Tq ~ cr_ ^ 0. The one-dimensional values 
(skyblue-circled dots) always lie on cr_. Note that in the 313 region lj > 0.53 for DVP, we can see 
''■o(^) ^ O(IO^); which implies that the contribution from the nonrenormalisable term in Eqs. ([28l IB9p . 
i.e. /S^ct" < 0(10-3) < 0(1), 0{\K\), is negligible compared to other terms in Eqs. (Hillll. Hence, 
our analytic solution still holds in the limit uj ^ 0(1) as discussed in See llll Cl 

Virialisation and characteristic slope: Fig. [S] shows the Q-ball properties plotted against the ra- 
tio ol S /U where S and U are the surface and potential energies (top panels), and the characteristic slope 
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FIG. 3: The top two panels show the three-dimensional numerical slopes —aja for two typical values of to for both 
DVP (left) and NDVP (right). The raw numerical data (red-solid and blue-dotted lines) matches continuously on 
to the analytical asymptotic profiles for large r (green-dashed and purple-dotted-dashed lines). The linear lines 
correspond to the Gaussian tails in Eq. psp where we can see the large shifts in the thin-wall limits of cj. The 
middle and bottom panels show, respectively, the hybrid profiles Eq. (|64|) and the energy density configurations 
for the various values of lo and D. The spikes of the energy density configurations exist in the DVP case but not 
in the NDVP case. 



Eq/loQ (bottom panels). For the DVP case where the thin-wall Q-ball satisfies ctq ^ cr+ it appears to 
be heading towards S /U ^ 1 &s u) ^ uj- = (sec Eq. p^ ). in all three cases. Also we predict that 
the thin- wall Q-ball in NDVP has S /U ^ Q (see Eq. (|14p ) and that it is consistent with what can be 
seen in the top right panel around cj = cj_ = 1. The bottom panels show analytically and numerically 
the characteristic slopes Eq/ojQ in both the thin and thick- wall limits. The analytic thin- wall lines 
(purple-dotted line for 2D and blue-dotted line for 31?) based on Eq. ((26)) are well fitted for the NDVP 
case with the corresponding numerical dots (red plus-dots for 2D and green crossed-dots for 31?) as long 
as (To — cr+i see the criteria in Fig.lH For the DVP case, our numerical data is seen to be heading in 
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FIG. 4: The initial value ao{uj) = <j{0) is plotted against uj. In the two panels the black-dashed and orange 
dotted-dashed lines show a±{uj), and these lines become closer for uj = ui- for both types of the potentials DVP 
(left, UJ- = 0) and NDVP (right, ui- — 1). Since ctq — cr+ for D = 2, 3 in the region uj ~ uj- where a+ ~ 1.28 in 
NDVP and a+ ~ 1.91 x 10^ in DVP, our analytical results in Sec lIII Bl are valid in this region. 



the right direction. The numerical solutions for both cases in the thick-wall region are well fitted by the 
analytic solution in general D given by the orange-dotted-dashed lines, in the second relation of Eq. ([33|) 
or Eq. (Hm. Fr om the virial relation Eq. (fT3| for Z) = 1 , we can only predict the extreme values of the 
ID characteristic slope, 7, in either the DVP or NDVP case once wc know what S/U is. To obtain that 
we rely on the numerical simulations and from the top two panels in Fig. [5l wc see that for the DVP 
case with D = 1, S/U appears to be heading towards unity, implying 7 » 1 in Eq. p^ . whereas for the 
NDVP case S/U ^ 1, implying 7 ^ 1 in Eq. ([T^ . Comparing these with the bottom two panels we see 
the behaviour for 7 appears to follow these predictions. 

Q-ball stability: Fig. [6]shows plots for both the classical (top panels) and absolute stability (bottom 
panels) with the stability threshold lines (black-dashed) for the cases of DVP (left) and NDVP (right). 
Let us consider the classical stability case first. For the thin- wall regime in DVP, notice that the numerical 
data of (rcd-dot-circlcs for 2D and grccn-dot-crosses for 3-D) are heading towards the analytic lines 
of Eq. ((TT]) . For the thick- wall case, on the other hand, the analytical lines of Eq. (orange-dotted- 
dashed) fit excellently with the numerical data in all dimensions, because Eq. (j34p is independent of D. 
Furthermore, the Q-ball is classically stable over all values of cu except for the ID thin- wall case where 
our analytical work cannot be applied. We saw this feature of unstable ID thin-wall Q-balls for the case 
of polynomial models in [s^]. For the absolute stability in the bottom panels, the analytical lines using 
Eq. and Eqs. IBlSp arc matched with the numerical dots for both the thin and thick-wall limits. 
Here, wc note how well the three-dimensional Q-ball (and also the higher dimensional ones as predicted in 
[37} ) can be described simply by our thin and thick- wall Q-balls. As our parameter set satisfies Eq. ([55)1 . 
we can see that absolutely stable Q-balls exist in DVP near the thick- wall limit. Because of the choice 
of UJ- = 1, the Q-ball in the NDVP case, however, is always absolutely unstable and most of the features 
are similar in terms of D. The analytical lines (top-right panel) in NDVP agree with the corresponding 
numerical data qualitatively better than the lines for DVP. 

To sum up our discussion of the Gravity mediated model, our analytical estimates of the characteristic 
slope and other properties of the Q-balls are well checked against the corresponding numerical results, 
even though we set a "fiatter" potential with \K\ = 0.1 < 0(1). 



B. Gauge-mediated potential 



This subsection presents numerical results showing the properties of gauge-mediated Q-balls with 
m = 1, = 2 in Eq. (|¥T|) . Although we have obtained analytical results for the potential, Eq. ([301), 
the potential is neither analytic nor smooth for all a. Therefore, we shall use the approximate potential, 
Eq. ((1T|) . see Fig. [2] and we expect that Eq. ([1T|) is a suitable approximation especially for the thin- wall 
limit UJ and large D. We will also see and explain the expected discrepancies that exist between the 
numerical and analytic results. 



17 



S/U 1.0 




1.0 



S/U 



0.6 



0.4 



0.0 



NDVP with K=-0.1 



ID 
2D 
3D 



0.2 0.3 0.4 0.5 0.6 



1.00 1.02 1.04 1.06 1.08 1.10 1.12 1.14 1.16 1.18 1.20 



3.0 



2.5 



i 2.0 



1.5 



ID 
2D 

3D 

thick wall analytic line 
2D analytical line 
3D analytical line 



DVP with K=-0.1 



*I?SSl:,»»^... 



0.2 0.3 0.4 0.5 0.6 

m 



1.14 
1.12 
1.10 
1.08 
1.06 
1.04 
1.02 
1.00 



II) 
21) 
31) 

thick wall analytic line 
2D analytical line 
3D analytical line 



NDVP with K=-0.1 



FIG. 5: The top panels show the ratio S/U where S and U are the surface and potential energies, and the bottom 
panels show the numerically obtained characteristic slope Eq/ujQ, in ID (skyblue circled-dots) , 2D (red plus-dots) 
and 3-D (green crossed-dots) . For comparison, in the bottom panels, the thin- wall analytic lines obtained using 
Eq. (|26p are also shown (purple-dotted line for 2D and blue-dotted line for 3-D) as are the thick-wall analytic lines 
obtained from Eqs. ([Ml [BTs]) ( orange-dotted-dashed for all D). The analytic lines match well with the numeric 
data in the appropriate limits, especially for the NDVP case. 



Hybrid profile: As we saw in earlier examples the numerical profiles we have obtained have errors 
for large ?', which correspond to cither undershooting or overshooting; thus, wc replace the numerical 
data in that regime by the exact asymptotic analytic solutions we obtained using the second relation of 
Eq. (j63p to smoothly continue the numerical solutions to the corresponding analytical ones. The hybrid 
profile in this model is 



'^W^ , ^/iW^(^-l)/%-rn„(.-i^„_) ^r i? <t<R ^^^^ 



where cr„„m is the numerical raw data, Rnum is determined by \^^^ +'mui + {'^num/'^num) I < 0.001, and 
we have again set Rmax = 60. We have calculated the following numerical properties using the above 
hybrid profile, Eq. ([55)) . up to -D = 3. 

Profile and energy density configuration: Fig. [7] shows the three-dimensional numerical slopes 
—a' I a for two values of uj (top), hybrid profiles (left-bottom) as in Eq. ([65]) . and the configurations 
for energy density (right-bottom). In the top panel, the raw numerical data (red-solid and blue-dotted 
lines) is matched smoothly onto the continuous asymptotic profiles Eq. (|65p for large r (green-dotted and 
purple-dashed lines). By fixing the numerical raw data using the technique Eq. ([65)) . wc show the profiles 
for various values of lo and -D, see the left-bottom panel. Also the peaks of the energy density cannot be 
observed in the whole range of see the right-bottom panel. 

Characteristic slope: In Fig. [51 we plot both the numeric and analytic characteristic slopes Eq/loQ 
(orange-dashed line for ID and blue-dotted line for 3-D). By substituting Eqs. ([5T1 [52)1 into Eqs. ([50l [Tl). 
wc have obtained the analytic slopes covering the whole range of lo. The 3-D analytic line agrees well with 
the numerical data except near the thick-wall limit. Similarly the 1-D analytic line agrees well only in the 
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FIG. 6: Classical stability for the top panels and absolute stability for the bottom panels for both DVP (left) and 
NDVP (right). The black-dashed lines indicate the stability thresholds for both classical and absolute stability 
in all panels. Q-balls found below the lines are stable either (both) classically or (and) absolutely. In the top 
panels, the analytical lines using Eqs. (|27l I34|) agree well quantitatively with the corresponding numerical data for 
thick-wall regimes, but not well in the thin-wall regimes. However the numerical plots look qualitatively similar 
to the analytical lines in the thin- wall limit as seen in the polynomial models [331. addition, the analytical 
lines for EqjmQ using Eqs. (|26l IB13]l match the numerical lines for both the thin and thick- wall limits. 



thin- wall limit. The origin of the discrepancies in the analytic versus numerical fits are the differences 
between the potentials themselves [Eq. (|40| in the analytical section (Sec. IIVP and Eq. (|4T|) ]. These 
differences are largest between 1 3 which in turn affects the region around 0.9 < w < 1.0, see 

Fig. [2] and Fig. [H 

Q-ball stability; Fig. [9] illustrates the stability of Q-balls: classical stability in the left panel and 
absolute stability in the right panel. The black-dashed lines in both panels indicate their respective 
stability thresholds where Q-balls under the lines are stable. We calculate the analytic lines for = 1, 3 
by substituting Eqs. (|51[ I52p into Eq. ((50)) and differentiating it with respect to io. The 31? numerical 
data can be matched with the analytic lines in both the thin and thick- wall limits. As in Eq. the 
three-dimensional Q-ball in the thick-wall limit is classically unstable. The numerical thick-wall Q-ball in 
IZ? is classically stable which differs from the prediction in Eq. ((6T|) . In the right panel, the analytic line 
for D = 3 agrees with the numerical data except in the thick-wall limit where the analytical lines for both 
\D and 3D do not match the corresponding numerical data. Furthermore, the thick-wall Q-ball in \D is 
absolutely unstable as predicted analytically in Eq. (|62p . but this fact cannot be observed numerically. 
The reasons for this discrepancy are as before a problem with our choice of potentials. We can see that 
the thin-wall Q-balls for any D are both classically and absolutely stable. 

To recap, our numerical results in the Gauge mediated case are generally well fitted by our analytical 
estimations. Observed discrepancies between the analytical predictions and numerical data arise from 
the artifact of our approximated smooth potential Eq. (|4T|) for the generalised Gauge mediated potential 
Eq. (PD|) . We have confirmed that the thin- wall Q-balls for any D are both absolutely and classically 
stable. 
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FIG. 7: The top panel shows the three-dimensional numerical slopes —a' /a for two values of uj. The raw numerical 
data (red-solid and blue-dotted lines) matches smoothly to the corresponding analytical asymptotic profiles for 
large r (green-dotted and purple-dashed lines). Both the left- and right-bottom panels show, respectively, the 
hybrid profiles Eq. (|65|l and the energy density configurations for the various values of u> and D. The spikes of 
energy density configurations do not exist even in the thin- wall limits. 



VI. CONCLUSION AND DISCUSSION 



We have explored stationary properties of Q-balls in two kinds of flat potentials, which are the gravity- 
mediated potential, Eq. (jisp. and the generalised gauge- mediated potential, Eq. ([40]) . Generally, the 
gauge-mediated potential is extremely flat compared to the gravity one; therefore, we cannot apply our 
thin-wall ansatz Eq. (|^ to the gauge-mediated case. By linearising the gauge-mediated potential, we 
obtained the analytical properties instead. For both potential types, we both analytically and numerically 
examined characteristic slopes as well as the stability of the Q-balls in the thin and thick-wall limits. Our 
main analytical results are summarised in Fig. 1101 

This present paper is of course related to our previous work (37j . The key differences arc that in the 
present work on thin-wall Q-balls we are assuming the value of cr+(tf) for the thin- wall limit uj ~ w_ 
depends weakly on lu and we have replaced the assumptio n a( Ro) < (J-^lo) by the equivalent assumption 
(made by Coleman) J7„ ~ J7^_ in the Q-ball shell region These in turn are related to the previous 
requirement that the surface tension t depends weakly on uj, which can be translated into the main 
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FIG. 8: The numeric characteristic slopes Eq/loQ and the analytic lines (orange-dashed line for ID and blue- 
dotted line for 3-D) which are calculated using Eqs. (|51l I52p in the whole range of ui. The 3D analytic line agrees 
with the numeric data well except near the thick-wall limit. Similarly the ID analytic line agree well only in the 
extreme thin- wall limit. 
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FIG. 9: The stability of Q-balls - Classical (left panel) and absolute (right panel). The black-dashed lines in the 
two panels indicate the stability thresholds for both classical and absolute stability where Q-balls under the lines 
are classically /absolutely stable. The analytic lines for D — 1, 3 are calculated by substituting Eqs. (|51l I52|l into 
Eq. (|5Up and differentiating it with respect to lo. 



assumptions: Rq 3> S, l/ii,aa ~ cr+, and ~ Ui^_ in the shell region. Furthermore, our analytic work 
agrees well with the numerical results for small curvature ^ with \K\ = 0.1; however, it is not clear that 
our analytic framework still holds even in the case of \K\ <^ 0{1), which corresponds to a case where the 
potential is extremely flat, see Eq. ((2T|) . 

Q-balls in gravity-mediated potentials: It is possible to obtain absolutely stable Q-matter with 
a small coupling constant. Kg. (|20p , for thcTlionrenormalisable term in Eq. ([TS]). For \K\ 0{1), a 
gravity-mediated potential cannot be really considered as flat, which allows us to apply our previous 
results, Eqs. ((26l [27|) . in [sj to describe the thin- wall Q-ball where ctq — <^+- In the thick- wall limit 
by reparameterising parameters in S^j and neglecting the nonrenormalisable term under the conditions 
S \K\ ^ have obtained the stationary properties of the Q-ball. We showed that the 



21 



thick-wall Q-ball is classically stable, and demonstrated that under certain conditions Eq. (|35p it can be 
absolutely stable. Although this analysis is much simpler than the analysis associated with imposing a 
Gaussian ansatz developed in Appendix [B] the former analysis assumed that the nonrenormalisable term 
is negligible at the beginning of the analysis. In the latter analysis, we have kept all terms in Eq. ([T7|) 
and shown that the nonrenormalisable term is indeed negligible in the limit to > C (to). Our results, 
Eqs. ((34l IB13|) , for the thick- wall Q-ball have recovered the previous results obtained in [s^, without 
any contradictions for classical stability conditions as opposed to the case of using a Gaussian ansatz in 
a general polynomial potential in which we showed that the ansatz led to a contradiction and corrected 
it by introducing a physically motivated ansatz [37| . This is because the Gaussian ansatz, Eq. (|B1[) . 
becomes the exact solution, Eq. (|Aip . in the Gravity mediated potential in the limit oj > 0{m) where the 
nonrenormalisable term is negligible. In Figs. [51 El the analytical lines agree well with the corresponding 
numerical plots in both the thin-wall and thick-wall limits. Under our numerical parameter sets, the 
Q-balls in DVP arc both classically and absolutely stable up to w < to, while all of the Q-balls in NDVP 
are absolutely unstable because of our choice, oj- = m. We believe that an absolutely stable Q-matter 
exists in NDVP when we take uj- < to. Since the Q-balls in both potential types are always classically 
stable, as can be seen in the top two panels of Fig. [HI except for the case of ID Q-balls in the thin-wall 
limit to which our analytical work cannot be applied since it holds only for > 2. We have also found 
the asymptotic profile Eq. for all possible values of ui, see the top two panels in Fig. [31 

Our analytical estimations on the value of do not agree well with the numerical results, because 
(To 7^ Nevertheless the other analytical properties are well fitted especially in NDVP, see bottom 
panels in Figs.|5l|6l The DVP in Eq. (fTS]) for small |iir| is extremely flat as the gauge-mediated potential 
in Eq. ([40]) where both of the potentials have uj- ~ 0. Notice that the asymptotic profile for the former 
case has a Gaussian tail, while the latter profile is determined by the usual quadratic mass term, see 
Eqs. ((38l [63l) . By assuming that the shell effects are much smaller than the core effects in the thin-wall 
limit, the difference of the tails can be negligible. Indeed, we can see the thin-wall numerical lines for 
both the classical stability and the characteristic slope look qualitatively and quantitatively similar to 
each other, as can be seen in both the top/bottom left panels of Fig. [Bland the panels of Fig. [9l Notice 
that the spikes of energy density in the gauge-mediated potential cannot be seen even though ~ 0, 
see Fig. [3 

Furthermore, wc know that the potential Ugrav can be approximated by ^TO^A/^l-'^lcr^"^'^! for small 
\K\ <C 0{1), then the potential in Eq. (jTSl looks similar to the confinement model in [isl . [l^. By 
neglecting the nonrenormalisable terms in the thick- wall limit, we can easily obtain the characteristic 
slope, 7 ~ 2+ i g i (ZJ-2) — by following the same technique as in Eq. (|TOl) . which does not depend 
on UJ but docs depend on D and It follows that Eq cx Q^^'' from Eq. (fTT|) . This result is obviously 
worse than our main results in Eqs. (|331 IB13p . see bottom two panels in Fig. [51 because we know that the 
Gaussian ansatz Eq. (|Bip can be the exact solution Eq. (jAl[) for U = Ugrav', thus, it is not so powerful 
to approximate Ugrav by ^rn^ M^^^^a^^^^^^ for small \K\. 

Q-balls in gauge-mediated potentials: For the Gauge mediated potential in Eq. (|40ll . we ob- 
tained the full analytic results in _D = 1, 3 over the whole range of io using Eqs. (jSTl [5^ . sec Figs. [SI [51 
In the thin-wall limit for niujR, (jjR » we reproduced the previously obtained results, Eq. (|54l) . 

in [U, m, 115 and showed that they are both classically and absolutely stable in Eqs. (|M1[551) . The one- 
and three-dimensional thick-wall Q-balls, on the other hand, are neither classically nor absolutely stable, 
see either Eqs. (|5T1 or Eqs. ([551 (BO)) , respectively. Since the potential, Eq. PU)) . is not differentiable 
everywhere, we have used the approximate potential, Eq. (|4T1) . instead in the numerical section. Figs. [3 [9l 
show that the numerical results agree with the analytical results in the thin- wall limit. The numerical 
data near the thick- wall limit and/or in the ID case differ from the analytic lines since the profiles 
are computed in the region where the two potentials between Eq. ([HI) and Eq. ((301) ai'e different, see 
Fig. [2l This differences come from the artifact of our approximated smooth potential Eq. ([4T1) against 
the generalised gauge-mediated potential Eq. ([3D]). 

The 3-D Q-balls: Although we have shown Q-ball results for an arbitrary number of spatial dimen- 
sions D, only three-dimensional cases are phenomenologically interesting. Q-balls in flat potentials give 
the proportional relation Eq oc Q^^^ , where 7 generally depends on D. The actual values of I/7 for 
three-dimensional thin- wall Q-balls are |, 1, and | in DVP, NDVP of gravity- mediated potentials and in 
gauge-mediated potentials respectively. It implies that the gauge-mediated Q-balls would be formed in 
the most energetically compact state for a large charge Q, so it is likely that such formed Q-balls would 
have survived any possible decay processes and thermal evaporation until the present day, and possibly 
become a dark matter candidate |29| . 
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Dynamics and cosmological applications: The dynamics of a pair of one-dimensional Q-balls 
has been recently analysed using momentum flux [sot - For a large separation between the Q-balls, the 
proflles develop the usual exponential tail, e"™'^'', in general polynomial potentials and in [s^l the authors 
showed that there was a solitonic force between them. Profiles in the gravity-mediated models and other 
confinement models, however, have different asymptotic tails, which may affect the detailed dynamics 
and the Q-ball formation [5l|, IH, [H, [13] . 

In a cosmological setting (thermal background), SUSY Q-balls are generally unstable via evaporation, 
diffusion, dissociation, and/or decay into todays baryons and lightest supersymmetric particles, if the AD 
field couples to the thermal plasma, which are decay products from inflaton, and/or if the field possesses 
a lepton number for the MSSM flat directions [U, [30|. Following our detailed analytical and numerical 
analyses of both gravity-mediated and gauge-mediated Q-balls, it is clear that this whole area of dynamics 
and cosmological implications of these Q-balls deserves further analyses. 
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APPENDIX A: AN EXACT SOLUTION 



In this appendix we will show that a Gaussian profile is an exact solution of the Q-ball equation in 
Eq. ^ with JJ^j = Ugrav — i^^o-'^ ill Eq. (fT5|) . Notice that the potential Ugrav becomes negative for 
gi/2|^f|^ < (J, hence the system is not bounded from below. The additional contribution from the 
nonrenormalisable term Unr compensates the negative term and supports the existence of Q-balls in the 
system. Although the Gaussian exact solution is no longer a solution for the full potential Ugrav + Unr 
in Eq. (|15|) . the solution we will obtain here provides hints in suggesting a reasonable ansatz for the 
thick- wall Q-ball as we will see later. 

Let us consider the following Gaussian profile: 

c^soi[r) = exp ( I , (Al) 

where we will see that m, M, and \K\ are the same parameters as in Eq. (|15p and p^j will be shortly 
determined in terms of the underlying parameters. By substituting Eq. (jAip into the left-hand side of 
Eq. ^ it leads to 

a 



Ugrav = —<j'{l^\K\\ni^—] ] (A2) 
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and 

where we set the integration constant as zero. Recall rn^ = rri^ — up' . Note that the constant M 
has the same mass dimension, (D — l)/2, as cr so that the only physical case is D = 3. The profile, 
Eq. (|A1[) . is an exact solution for Ugrav with the "core" radius Rq = ^2lm?-\K\ [s^, which is very large 
compared with for small \K\ ^ and satisfies the boundary conditions for Q-balls, namely 

(t'(0) = = (t(oo) = u'ipo) [37|- In the extreme limit a; ^ m, we obtain ^ for \K\ < 0(1) 
which implies ctq = o'(O) — s- 0. For large cr, the potential becomes asymptotically flat, tending towards 
an infinite negative value. By adding the nonrenormalisable term Unr, the potential Ugrav is lifted for 
large a in Eq. (fT5|) . then the full potential Ugrav + Umr is bounded from below, see Sec. IIII Al Wc can 
see the ansatz given in [s^ corresponds to the case where p^j — M, which is valid only for \K\ ^ 0(1) 
and u ~ m, see Eq. (|A3|) . 



APPENDIX B: THICK- WALL Q-BALL WITH A GAUSSIAN ANSATZ 

In this appendix, wc will investigate the thick-wall Q-ball in gravity-mediated models by introducing a 
Gaussian ansatz and keeping all terms in Eq. pT]) as opposed to the analysis in Sec lIII Cl By using this 
profile we can perform the Gaussian integrations, and will obtain the generalised results of Eqs. ([33l [34]) 
in Sec. IIII CI The test profile for the case, w > 0(m), coincides with the solution a sol in Eq- (|Aip . which 
implies that the nonrenormalisable term U nr in Eq. psp is negligible. 

To recap, the notation we have adopted in Eq. pT|) is ct = cr/Af, ui = uj/m, (3'^ is defined in Eq. p6)) 
and we are considering the case of n > 2. To begin with we introduce a Gaussian ansatz inspired by 
Eq. ([AT]) for the potential Eq. 

a(r) - A^cxp(-k2^V2), (B1) 

where ctq = o'(O) = A^^^ finite, and A^^, will be functions of uj implicitly. A^^^ should not be confused 
with the coupling constant A in Eq. (jl5[) . Both A^^; and can be determined by cxtrcmising the Euclidean 
action S^] hence the actual free parameter here will be only to. It is crucial to note that A^^ cannot be 
infinite in the thick-wall limit since wc know that A^j is finite and tending to 0. If the nonrenormalisable 
term Unr^ is negligible, we can expect A;^ ^ Puj/M ~ cr_(a;) and ~ |i^|m^ due to Eq. (|A1[) . which 
implies that the "core" radius Rq of the thick- wall Q-ball is Rq ~ y^2/m'^\K\. For the extreme thick- wall 
limit a; 3> TO, wc shall also confirm A(j — > 0, which means (Tq — > 0. 

By substituting Eq. (jBip into Eq. ([T|) with the potential Eq. (fTT)) . we obtain Q and S^j using the 

following Gaussian integrations: flD-i drr^^^e^'^^ = (f)^^^ for real k where flo-i = r{D/2) • 

Q = Af^TT^/^^A^^J^, (B2) 

= Af^TT^/^/^-^ [A(a«^, X^) + B{cj, A„) + C(A^)], (B3) 
r)\2 

where Ain^, \^) = —^{til + \K\m^), (B4) 



Thus, 



4 

Bico, A^) = ( 1 - ^ - m\ InA^ ) , C(A^) = m'p'XZ ( - ) . (B5) 



TO^A^ / iu^ \ . /2^^/' 



n 

Notice that A(k^, X^^) comes from the gradient term and the logarithmic term in S^, and depends on 
both and A^^. Similarly, B(uj, X^) is given by the quadratic term in the potential Eq. (|17p and depends 
both on Atj and explicitly on uj, whereas C{X^) arises simply from the nonrenormalisable term in the 
potential. An alternative (but in this case more complicated) approach to obtain Q would be the use of 
Legcndrc transformations in Eq. pOp . 

By cxtrcmising in terms of the two free parameters k^^j and A^^ : 

1^ = 0, ^=0, (B6) 

we obtain 

A + B + C^^Y"^, A + B + — ^ (B7) 
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which impUes that 



= \K\ - {n - 2)(3^X: 



D/2 



> 0, 



(B8) 



where we have ehminated the A + B terms in the two expressions of Eq. (jB7[) . Using Eq. (jB8[) and the 
second expression of Eq. (jB7[) , wc obtain the relations between to and 



(B9) 
(BIO) 
(Bll) 



m 



duj 



l + |X|(i:'-l-2hiA„) for ~ 0(1), 
l-2|is:|lnA^ for liCl < C(l), 

< 0, 



|X|m2F IKW 



where we have differentiated Eq. (jB9p with respect to oj to obtain Eq. (jBlip and have defined F as 



F =1- {n 
that both K, 



2) 



2{n+D)-nD jl'^ 



A"; 



f 



2(n+D)-nD 



\K\). Equations (|B8l [B9| imply 



and Ai^ are functions of oj; however, these arc not solvable in closed forms unless the particular 
limits, which were introduced in Sec lIII Cl are taken, as we will now show. Comparing Eqs. (|B9|, IB11|) with 
Eqs. ([28l [29 ll . wc can see an extra contribution of C(|iir|) in Eq. (|B9[) . which is not present in Eq. ([28|) . 
This difference of {D — \)\K\ arises because in calculating Eq. (jB9[) we have used A^^, whereas we have 
used cr-{ijj) in obtaining Eq. ((28|) . and although related they are not precisely the same. In the extreme 
thick-wall limit ^ m, and from Eq. (jB9P this implies A^^ 0+ (recall from Eq. (|Bip that A^^ has to 
remain finite). Considering the nonrenormalisable term in Eq. (|B9p . the fact that 0^ < |/\ | < 0(1) and 
A^^ 0+ with n > 2. implies that this term is subdominant and can be ignored. As long as A^^ < 0(1), 
then F ^ \ and the second relation of Eq. (jBlip follows, which implies that A^^; is a monotically decreasing 
function in terms of lo. The limit A^^ ^ 0(1) corresponds to > 0{m), see Eq. (|B9p . We will call this 
the "moderate limit" and represent it by The other case, oj ^ m (or equivalently ^ 0(1)), we 
shall call the "extreme limit" and represent it by Depending on the logarithmic strength of \K\, we 
can obtain Eq. (jB10[) . which leads to the approximated expressions for A^^ and can also obtain from 
Eq. (HH) 



(T^{uj) for |A| < 0(1) m'^ 



(B12) 



where is independent of lu in both the "moderate" and "extreme" limits. 

Using Eqs. (|B2| IB3[1 and Ec^. (jB7p . we obtain the characteristic slope in both the "moderate" and 
"extreme" limits. 



Eq 

LUQ 



= 1 



2^2 



nf\K\ 
2^2 



1. 



(B13) 



In order to show their classical stability, wc shall differentiate Q with respect to to using Eqs. (|B8I IB9 
and Eq. (|BTT|) : 



LO dQ 
Q duj 



dtu \ Q 



m'^\K\F 

'm?\K\ 
1 

2^ 



m 



\K\) 



2u^ 



2^2 



'm?\K 
(n- 2)cj2 
^ 'm?\K\F 

1 > 0, 



<o. 



\K\) 



(B14) 



(B15) 



where we have taken the "moderate limit" and "extreme limit" and used ^ to^IXI, F = 1 



2{n+D)-nD / 2 



I if I) ~ 1. The classical stability condition Eq. (|B14[) is consistent with Eq. (|B15[) . 
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and is consistent with Eq. p^ . This is different from the result we obtained for the polynomial poten- 
tials [see Eq. (74) in [37|], because in that case the Gaussian ansatz does not give the exact solution 
unlike here in Eq. (jBip where it does become the exact solution Eq. (|Aip in both limits. The results, 
Eqs. (jB13[ [BT4|) and Eq. (jB15p . in both the "moderate" and "extreme" limits recover the key results, 
Eqs. ([551 IM)) . and arc independent of D; hence, the thick-wall Q-balls for all D have similar properties. 
We can also see the small additional effects arising from the nonrcnormalisable term in Eqs. (jB141 IBT5|) . 

Let us summarise the important results wc found in this appendix. By introducing a Gaussian test 
profile Eq. (|Bip inspired by the exact solution Eq. (jAip for Ugrav, we computed the Euclidean action 
and the charge Q using Gaussian integrations. Then, we extremised in terms of A^^ and in 
Eq. (|B6[) . which gave the relations of both X^^ and as a function of lu. By introducing two limits called 
"moderate limit" and "extreme limit", we confirmed that the ansatz, Eq. (jBip . approaches Eq. (jAip 
in the "moderate limit" . We established that the results Eqs. (jB131 IB14I) and Eq. (jB15[) recovered the 
previous results in Eqs. (j33i I34[) which are obtained simply by reparanicterising in S^^ and extracting the 
explicit oj-dependence from the integral in Suj with U = Ugrav where the nonrcnormalisable term was 
neglected at the beginning of the analysis by applying L'Hopital rules. 

In addition, we would like to emphasise the main differences between our work and other earlier analyses 
in the literature [lO,!!^- The analytical framework adopted in [i^ is valid only for \K\ = 1, D = 3, n ~ 4. 
Our work has shown that this can be generalised to arbitrary integer values of D and n(> 2) under the 
conditions < \K\ < 0(1), and that the thick-wall Q-ball can be classically stable. In Scc lIII Cl we also 
found that the thick- wall Q-ball may be absolutely stable under certain additional conditions, Eq. (pS)) . 
Furthermore, Enqvist and McDonald in [s^] analytically obtained the same "core" size of thick-wall Q- 
balls, although they obtained a slightly different value for Eq/Q (see their Eq. (112)). The reason for 
this is because their ansatz assumed A^^ ~ 1 in Eq. (|Bip by simply neglecting the nonrcnormalisable 
term, which implies that the third term of B{ijj, A^^) and term C(A(^) in Eq. (|B5[) are absent. Hence, 
their analysis is vahd for \K\ <^ C(l) and uj ~m, see Eq. (|A3|) . We, however, have kept all the terms in 
Eq. P7|) and used a more general ansatz, which can be applied for \K\ < 0(1) and w > 0{m) with the 
restricted coupling constant of the nonrcnormalisable term /3'^ ^ \K\. In summary, in this appendix we 
have extensively investigated analytically both the absolute and classical stability of Q-balls in Eq. (jB13p 
and Eq. (|Bl4)) . 
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